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The magnetically ordered, low temperature phase of Ising ferro- magnets is manifested within 
, the associated Fortuin-Kasteleyn (FK) random cluster representation by the occurrence of a single 

positive density percolating cluster. In this paper, we review our recent work on the percolation 
signature for Ising spin glass ordering — both in the short-range Edwards-Anderson (EA) and 
I infinite-range Sherrington-Kirkpatrick (SK) models — within a two-replica FK representation and 

also in the different Chayes-Machta-Redner two-replica graphical representation. Numerical studies 
of the ±J EA model in dimension three and rigorous results for the SK model are consistent in 
i-^ ' supporting the conclusion that the signature of spin-glass order in these models is the existence of a 

single percolating cluster of maximal density normally coexisting with a second percolating cluster 
of lower density. 
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I. INTRODUCTION 



> 

The question of whether laboratory spin glasses — or the theoretical models used to represent them — have a 
I thermodynamic phase transition remains unresolved despite decades of work [11, 0! 3 ■ Although the infinite-range 
' ^ Sherrington-Kirkpatrick (SK) Ising spin glass is easily shown to possess a phase transition [1, Q, the existence of 
Q ■ one in the corresponding short-range Edwards-Anderson (EA) [6] Ising model (on the cubic lattice Z'^) has not been 
^ [ established in any finite dimension. While some evidence for a transition has been uncovered through high-temperature 
'— exp ansions 0,@|, analytical studies of variable long-range Id models [ll, and extensive numerical simulations (T], 
^ ' [ll[ [l^ , the issue remains unresolved [T§| . 
. Random graph methods, and in particular the Fortuin-Kastelyn (FK) random cluster (RC) representation p^. [l5|. 

■ provide a set of useful tools for studying phase transitions (more specifically, the presence of multiple Gibbs states 
^\ , arising from broken spin rotational symmetry) in discrete spin models. In these representations spin correlation 
t functions can be expressed through the geometrical properties of associated random graphs. FK and related models 

\ are probably best known in the physics literature for providing the basis for powerful Monte Carlo methods for 

■ studying phase transitions [l6l . IITI . [isj , but they have also proved important in obtaining rigorous results on phase 
CD ] transitions in discrete-spin ferromagnetic (including inhomogeneous and randomly diluted) models (e.g., [l^ [20j). 
OO Because of complications due to frustration, however, graphical representations have so far played a less important 

role in the study of spin glasses. 
J> [ The goal of the studies presented here is to construct a viable approach which uses random graph methods to address 
the problem of phase transitions and broken symmetry in spin glass models. In earlier papers j2ll. [2^. we studied 
the "percolation signature" of spin glass orderin g w ithin two different graphical representations — the two-replica 
model of Chayes, Machta and Redner (CMR) [23l. |23| and a two-replica version of the FK representation (TRFK), as 
proposed in Sec. 4.1 of [25)]. The purpose of that analysis was to show that FK methods could be utilized to study 
spin glass phase transitions. The result of this work was the uncovering of strong evidence that the existence of a 
spin glass transition coincides with the emergence of doubly percolating clusters of unequal densities. This scenario is 
more complex than what occurs in ferromagnetic models, where the phase transition coincides with percolation of a 
single FK cluster. 
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In what follows, we first review the FK random cluster representation for ferromagnetic models, and then describe 
our analytical and numerical results for spin glasses. 

II. THE FORTUIN-KASTELEYN RANDOM CLUSTER REPRESENTATION 

In this section we briefly review the Fortuin-Kasteleyn random cluster representation fT3, Ts^ which relates the 
statistical mechanics of Ising (or Potts) models to a dependent percolation problem. Our focus will be on Ising 
models throughout, but the analysis is easily extended to more general Potts models. 

A. Ferromagnetic Models 

We start by considering a nearest-neighbor Ising ferromagnet, whose couplings Jxy > are not necesarily identical: 

where, as already noted, Jxy > and the sum is over nearest neighbor pairs of sites in Z"^. Each such coupling is 
associated with an edge, or bond, (x, y), with the set of all such bonds denoted E''. 
The RC approach introduces parameters Pxy G [0, 1) through the formula: 

Vxy ^ I - eM-Wxy\] , (2) 

where /3 is the inverse temperature. One can then define a probability measure fiRC on {0, l}^'', that is, on 0- or 
1-valued bond occupation variables ujxy It is one of two marginal distributions (the other being the ordinary Gibbs 
distribution) of a joint distribution on O = { — 1,+!}^ x {0,1}^ of the spins and bonds together (such a joint 
distribution will be introduced in Sec. Ill Cjl . The marginal distribution /iRc is given formally by 

mciWxy}) = 2#«'^-» flindiWxy}) luiWxy}) , (3) 

where Zrc is a normalization constant, #({wxy}) is the number of clusters determined by the realization {ujxy}, 
fJ-indii'-^xy}) is the Bernoulli product measure corresponding to independent occupation variables with fJ.ind{{'-^xy — 
1}) = Vxy, and lu is the indicator function on the event U in {0, 1}^ that there exists a choice of the spins {<7x} so 
that Jxy(^xyO-xCry > for all {x,y) [H, [13, [111 . 

There are several things to be noted about ([3]). The factor 2 in the term 2'^(^"="«^) arises because we've confined 
ourselves to Ising models, so that every connected cluster of spins (each such connected cluster consists of all satisified 
bonds) can be in one of two states (in the ferromagnetic case, all up or all down); in a g-state Potts model, this 
term would then be replaced by q#(i'^^v}\ More importantly, the indicator function on U, which is the event that 
there is no frustration in the occupied bond configuration, is always one for the ferromagnet; consequently, this term 
is superfluous for ferromagnetic models. We include it, however, because it will be needed when we generalize to 
models with frustration. Finally, we note that finite-volume versions of the above formulas, with specified boundary 
conditions, can be similarly constructed. 

In the case of a ferromagnet, there are general theorems [29j which ensure that when percolation occurs, there 
is a unique percolating cluster. It then easily follows that RC percolation within the FK representation (or "FK 
percolation" for short) corresponds to the presence of multiple Gibbs states (in the ferromagnet, magnetization up 
and magnetization down) , and moreover that the onset of percolation occurs at the ferromagnetic critical temperature. 
To prove this, it is sufficient to show that FK percolation is both necessary and sufficient for the breaking of global 
spin flip symmetry in the ferromagnet. To see that FK percolation is a necessary condition, note that the contribution 
to the expectation of ctq from any finite RC cluster is zero: if a spin configuration a is consistent with a given RC 
bond realization within such a cluster, so is —a, and both will be equally likely. As a consequence, (cro) = in infinite 
volume in the absence of RC percolation. 

To see that RC percolation is a sufficient condition for the magnetization order parameter to be nonzero, consider a 
finite volume with fixed boundary conditions, i.e., a specification ct^; = ±1 for each ct^; G OAl- For the ferromagnet, 
by first choosing all Wx = +1 and then all Wx — —1, one can change the sign of the spin ctq at the origin even as 
L oo. That is, boundary conditions infinitely far away affect ctq: which is a signature of the existence of multiple 
Gibbs states. 
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When attention is confined to ferromagnetic models, the mapping of the FK formahsm to interesting statistical 
mechanical quantities is straightforward (and intuitive); for example 

(ax<Jy) = mc(a; ^ y) , (4) 

where (axCy) is the usual Gibbs two-point correlation function and fJ.Rc{x <-> y) is the RC probability that x and 
y are in the same cluster. Similarly, with "wired" boundary conditions (i.e., each boundary spin is connected to its 
nieghbors), one has 

{ax)+ = URcix ^ oo) . (5) 

So for ferromagnets, a phase transition from a unique (paramagnetic) phase at low (3 to multiple infinite- volume Gibbs 
states at large p is equivalent to a percolation phase transition for the corresponding RC measure. 



B. Spin Glass Models; The TRFK Representation 

For spin glasses (or other nonfcrromagnets with frustration) the situation is more complicated. Now for two sites 
X and y, ([4]) becomes 

{cTxCTy) ^ {lx^yV{x,y))Rc; v{x,y)= Yl sgn{Jx'y'), (6) 

{x'y')eC 

where C is any path of occupied bonds from x to y. By the definition of U, any two such paths C and C in the same 
cluster will satisfy U{x'y')eC sgn{Jx'y') = Il{x'y')eC' sgn(J:.'y')- 

Just as in the ferromagnet, if percolation of a random cluster occurs in the FK representation, the percolating 
cluster is unique (in each realization of FK spins and bonds) psl [30l |. In spite of this, RC percolation alone is no 
longer sufficient to prove broken spin-flip symmetry in the EA spin glass. This is because even in the presence of RC 
percolation, it is unclear whether there exist any two sets of boundary conditions that can alter the state of the spin 
at the origin from arbitrarily far away. Although the infinite cluster in any one RC realization is unique, different RC 
realizations can have different paths from ^ SA^, and because of frustration this can lead to different signs for ctq. 
So percolation might still allow for (cto) — > as L ^ oo, independently of boundary condition. 

Indeed, it is known that FK bonds percolate well above the spin glass transition temperature. For the three- 
dimensional Ising spin glass on the cubic lattice, Fortuin-Kasteleyn bonds percolate at /3fk,p ~ 0.26 [HI while the 
inverse critical temperature is believed to be /3c = 0.89 ± 0.03 [s!]. Near the spin glass critical temperature, the giant 
FK cluster includes most of the sites of the system. For this reason, the Swendsen-Wang algorithm, though valid, is 
inefficient for simulating spin glasses. 

However, it is an open question as to whether the presence of "single " (see below) FK percolation would lead to a 
slower, e.g. power-law, decay of correlation functions even though the Gibbs state is unique. If this were to happen, 
then the onset of single FK percolation would imply a phase transition in the spin glass, though not multiple Gibbs 
states and hence no broken spin- flip symmetry. This possibility was suggested in [25]; but to date, no evidence exists 
to support it. 

Nonetheless, single FK percolation remains a necessary condition for multiple (symmetry-broken) Gibbs phases in 
the spin glass, for the same reason as for the ferromagnet. A slightly stronger version of that argument [28[ proves 
that the transition temperature for an EA spin glass, if it exists, is bounded from above by the transition temperature 
in the corresponding (disordered) ferromagnet. 

The essential difference (from the point of view of FK percolation) between the ferromagnet and the spin glass is 
the factor Ijj in ([3]). It is somewhat easier to discuss this factor in the context of finite volumes. Let A^ denote the 
L'^ cube centered at the origin, and let A^ denote the set of bonds b = {x,y) with both x and y in A^. For models 
containing frustration, U is typically not all of {0, 1}"^^, unlike for the ferromagnet; it is the set of all unfrustrated bond 
configurations (so that the product of couplings around any closed loop in such a configuration is always positive). 

So is there a way to extend the above considerations to arrive at a sufficient condition for multiple Gibbs states in 
spin glass models? We begin by noting that FK clusters identify magnetization correlations; but spin glass ordering 
is manifested by the Edwards- Anderson (EA) order parameter, and not the magnetization, becoming nonzero. The 
EA order parameter qea can be defined with respect to two independent replicas of the system, each with the same 
couplings {Jxy}- Denoting the spins in the two replicas by {cx} and {tx}, each taking values ±1, qEA is defined in 
terms of the overlap. 



Q = N-^^GxTx , 
{x} 



(7) 
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in the limit as the number of sites N ^ oo. In general, Q is a random variable whose maximum possible value is Qea, 
but in the case where (in the limit N oo) {cr^} and {t^} are drawn from a single pure state, Q takes on only the 
single value qea- 

Using this as a guide, it appears that one possibility, proposed in [2ll. [2^. [25j . for extending FK methods to spin 
glasses is to use what might be called double FK percolation. Here one expands the sample space H, to include 
two independent copies of the bond occupation variables (for a given configuration), and defines the variable 
r^y = LUxyUj'^y^ where uj^y and uj'^y are taken from the two copies. One then replaces percolation of {i-^xy} in the single 
RC case with percolation of {rxy}- It is not hard to see that this would be a sufficient condition for the existence of 
multiple Gibbs phases (and consequently, for a phase transition). 

This is not the only way, however, to use "double" percolation of RC clusters in some form to arrive at a condition 
for spin glass ordering; the above describes what we denoted in the Introduction as the TRFK approach. In the 
next subsection we present a closely related approach, the CRM two-replica graphical representation, which requires 
a more lengthy description. Once this is done, we will discuss how both representations can be used to illuminate the 
nature of phase transitions in short-range spin glass models. 

C. The CMR Representation 

In describing these "double FK" representations, we will find it useful to reformulate our description of the relevant 
measures slightly using a joint spin-bond distribution introduced by Edwards and Sokal 33]. The statistical weight 
W for the Edwards-Sokal distribution on a finite lattice is 

Wi{ax}, Wxyhp) - pl'^l(l - p)^^-l"l A({a4, {ioxy}) . (8) 

Here = ^^(k j^) is the number of occupied bonds and Ni, is the total number of bonds on the lattice. The factor 
A({cra;}, {uJxy}) is introduced to require that every occupied bond is satisfied; it is defined by 

A({a.}, M) = { ; l^^^^Z^ ^y-- ^-^^'^^ ^ « (9) 

With p defined as in ^ and with all J^y = J > 0, it is easy to verify that the spin and bond marginals of the 
Edwards-Sokal distribution are the ferromagnetic Ising model with coupling strength J and the Fortuin-Kasteleyn 
random cluster model (cf. ([3])), respectively. 

We can now adapt the above representation to the ± J Ising spin glass. (With minor modifications, it can also be 
adapted to Gaussian and other distributions for the couplings.) The corresponding Edwards-Sokal weight is the same 
as that given in Eq. ([5]). The A factor must still enforce the rule that all occupied bonds are satisfied, 

A({aJ, {.xy}; {Jxy}) = { ; ''JZ"^^' ^-"-"^"^ - " (10) 

SO (cf. ([3])) Ijj — maxj-g-a,} A({(7a;}, {uJxy}] {Jxy})- The spin marginal of the corresponding Edwards-Sokal distribution 
is now the Ising spin glass with couplings {Jxy}- But the relationship between spin-spin correlations and bond 
connectivity is complicated by the presence of antiferromagnetic bonds. Now one has 

{ax(Ty) = (11) 

Prob{x and y connected by even number of antiferromagnetic bonds} 
— Prob{x and y connected by odd number of antiferromagnetic bonds}. 

The two- replica CMR graphical representation, introduced in [23,[23|j incorporates, in addition to the spin variables 
{cJx} and {tx}^ two sets of bond variables {oJxy} and {rjxy}^ each taking values in {0, 1}. Now the Edwards-Sokal 
weight is 

y^{{(^x}, {Tx}, {uJxy}-, {llxy}; {J'xy}) 

= Bh\nc{{tjJxy})Brcd{{rixy})^{{(^x}, {Tx} , {t^xy}\ {Jxy})Y{{<7x}, {Tx} , {Vxy}) (12) 

where the S's are Bernoulli factors for the two types of bonds, 

B^yU{^xy}) = - ^biuc)'^^-!'^! (13) 

i?rcd({^.,})=P|;l(l-^red)^'-|''l (14) 
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and the bond occupation probabilities are 



Pbiuc = l-exp(-4/3|J|) 
T'.ed = l-exp(-2/3|J|). 



(15) 
(16) 



So Bi^\uci{'-Oxy}) = fJ-indii'-^xy})- Thc A and r factors constrain where the two types of occupied bonds are allowed, 



Mi^x}, {Tx}, {^^xy}; {Jxy}) 

_ r 1 if for every xy: JxyUJxyC^xCry > and JxyCOxyTxTy > 
~ [ otherwise , 

^iWx}ATx}AVxy}) 

1 if for every xy: rj^yCr^ayT^Ty < 
otherwise . 



(17) 



(18) 



We refer to the tj-occupied bonds as "blue" and the jy-occupicd bonds as "red" . The A constraint says that blue 
bonds are allowed only if the bond is satisfied in both replicas. The F constraint says that red bonds are allowed only 
if the bond is satisfied in exactly one replica. 

It is straightforward to verify that the spin marginal of the two-replica Edwards-Sokal weight is that for two 
independent Ising spin glasses with the same couplings, 



yVi{ax}, {Tx}, {cJxy}, {Vxy}; Uxy}) 



const X exp 



P ^ ^ Jxyi^^x^y ~^ '^x^y 

{xy) 



(19) 



Connectivity by occupied bonds in the two-replica representation is related to correlations of the local spin glass 
order parameter. 



^x^x 



It is straightforward to verify that 



{QxQy) = Prob{x and y connected by even number of red bonds} 
— Probjo; and y connected by odd number of red bonds}. 



(20) 



(21) 



As in the case of the FK representation, a minus sign complicates the relationship between correlations and connectivity 
but in a conceptually different way. The second term in Eq. (|21|) is independent of the underlying couplings in the 
model and is present for both spin glasses and ferromagnetic models. 

As noted earlier, for ferromagnets (in the absence of non-translation-invariant boundary conditions), the signature 
of ordering is a single percolating FK cluster. For spin glasses, the situation is more complicated as there can be 
more than one percolating cluster. However, if the CMR graphical representation displays a single percolating blue 
cluster of largest density, one can similarly show broken symmetry, for either EA or SK models. This is because 
one can impose "agree" or "disagree" boundary conditions between those ax and boundary spins belonging to 
the maximum density blue network. In the infinite volume limit, these two boundary conditions give different Gibbs 
states for the a-spin system (for fixed r) related to each other by a global spin flip (of a). 

It is a separate matter, however, to relate in a simple way the EA order parameter to the density difference in blue 
(or more generally, doubly percolating) clusters. Intuitively, it seems that there should be a simple correspondence, 
and in fact such a relation is easy to show for the SK model: here the overlap Q is exactly the difference in density 
between two percolating blue clusters [2l| . But it is not immediately obvious that a similar density difference can be 
simply related to the EA order parameter (although the arguments above imply that here also a density difference 
implies a nonzero Q.) This is because in a two-replica situation, it is not immediately obvious that there is no 
contribution from finite clusters in short-range models. 

For the TRFK representation, similar reasoning shows that the occurrence of exactly two doubly-occupied perco- 
lating FK clusters with different densities implies broken symmetry for the spin system [25[ and that Q should equal 
(and once again, does equal in the SK model) the density difference. In \2l\ we presented preliminary numerical 
evidence that there is such a nonzero density difference below the spin glass transition temperature for the d = 3 EA 
±J spin glass (for both the TRFK and CMR representations). We will discuss these results further in Sec. IIVI But 
now we turn to a review of some rigorous results, particularly for the SK model, that appeared in and |22j |. 
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III. RIGOROUS RESULTS 

The SK model permits a fairly extensive rigorous analysis of both the CMR and TRFK representations which, 
when combined with other known results, permits a sharp picture to emerge of the connection between double FK 
percolation (we use this more general term to refer to any representation that relies on two FK replicas, such as 
CMR or TRFK), and a phase transition to a low-temperature spin glass phase with broken spin- flip symmetry. The 
surprising feature to emerge from this analysis is that in the CMR representation, there already exist well above the 
transition temperature two percolating networks of blue bonds, of equal density and in all respects macroscopically 
indistinguishable. Below the critical temperature Tc, the indistinguishability is lifted: the two infinite clusters assume 
different densities. 

On the other hand, there is no double FK percolation at all above Tc in the TRFK representation, but below Tc 
there are again two percolating double clusters of unequal density. 

We do not yet know whether this difference in the two pictures above Tc persists in the EA model, but numerical 
evidence (to be discussed in the next section) so far appears to indicate that it does not: in the EA spin glass, both 
representations seem to behave similarly to CMR in the SK model. 

The SK Hamiltonian for an A'^-spin system is 

Hn = — y= ^ J^j<T^aj , (22) 

l<i<j<Af 

where i,j = I, . . .N are vertices on a complete graph. The couplings Jij are i.i.d. random variables chosen from a 
probability measure p satisfying the following properties: 

1) The distribution is symmetric: p{u) — p{~u). 

2) p has no 5-function at w = 0. 

3) The moment generating function is finite: dp{u) e"*'"' < co for all real t. 

4) The second moment of p is finite; we normalize it to one: dp{u) v? — 1. 

So the results of this section hold for Gaussian and many other distributions (although not diluted ones, because of 
requirement (2)), but the analysis is simplest for the ± J model, to which we confine our attention for the remainder of 
this section. We therefore assume that Jy = ±1 with equal probability; for this distribution (or any other satisfying 
property (3)) /3e = 1 fli- ^ 

It is relatively easy to see, using a heuristic argument presented in [2l|, at what temperature a single FK cluster 
will form. The SK energy per spin, u, is given above the critical temperature by u = —[3/2. Therefore, for large N 
the fraction fs of satisfied edges is 

fs-\~ uN-''^ . (23) 

From ([2|) it follows that a fraction "Pfk = 1 ^ exp(— 2/3A^^^/^) « 2/3 A^^^/^ of satisfied edges are occupied. According 
to the theory of random graphs (see ^^), a giant cluster forms in a random graph of N vertices when a fraction x/N 
of edges is occupied with x > 1, and there is then a single giant cluster [s^. So if edges are satisfied independently 
(which of course they're not — this is why this argument is heuristic only), then single replica FK giant clusters 
should form when j3 = xN^^^"^ when x > 1. The single replica FK percolation threshold is therefore at 

/3fk,p - Af-i/2 ^ (24) 

and above this threshold, there should be a sing le g iant FK cluster. 

This simple argument can be made rigorous [25|, but the basic ideas are already displayed above. The rigorous 
argument obtains upper and lower bounds for the conditional probability that an edge {xoyo} is satisfied, given the 
satisfaction status of all the other edges. If these bounds are close to each other (for large A^) then treating the 
satisfied edges as though chosen independently can be justified. We now describe that argument. 

To avoid the problem of non-independence of satisfied edges, we fix all couplings Jxy but one, which we denote 
Jxoyoi ^^^d ask for the conditional probability of its sign given the configuration of all other couplings and all spins 
(7x- The ratio Z+/Z^ of the partition functions with Jxoyo = +N^^^^ and Jxoyo ^ ~N~^/^ satisfies 

exp(-2/3A^-i/2) < \Z+/Z_\ < exp{2f3N-'^^^) . (25) 
The conditional probabilities P± that Jxoyo ~ ±A^^^^^ therefore satisfy 

^-4p/VN < e-^/'/^lZ^/Z+l < P+/P_ < e^^^^\Z^/Z+\ < e^^l^ . (26) 
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Let Ps (Pu) be the conditional probability for any edge {xoj/o} to be satisfied (unsatisfied) given the satisfaction 
status of all other edges. These must then be bounded as follows: 



and therefore 



-4P/VN ^pjp^^^4p/VN ^ (27) 



^ - 0{/3/Vn) = {e^P/^ + l)-i < < (e-4/3/v^ + i)-i = 1 + 0[(i/VN) . (28) 



One now obtains rigorously the same conclusions as before — i.e., ([24]) is valid with a single giant FK cluster for 
(3 = (3n> xN~'^/'^ with any x > I. 

Essentially the same argument is used for the more interesting case of double percolation. Here there are two 
spin replicas denoted by a and r. We are now interested in percolation of doubly satisfied edges; the spins at 
the vertices of each such edge must satisfy UxTx = o'yTy (and then will be satisifed for exactly one of the two 
signs of Jxy)- We note an immediate difference between the CMR and TRFK representations: for the former, 
doubly satisfied edges are occupied with probability Vcmr = 1 ~ exp(— 4/3iV^^/^) ^ 4/3 A^^^/-^, while for the latter, 
^TRFK = [1 ~ exp(— 2/37V~^/^)]^ ~ 4/3^/iV. It seems likely that this difference occurs only for the SK model; the 
various factors of N are absent in the EA model. 

We can proceed much as in the single-replica case by dividing all (a, r) configurations into two sectors — the agree 
(where ax = Tx) and the disagree sectors (where ax = —Tx)- We also denote by Na and Nd the numbers of sites in 
the sectors and denote by Da — Na/N and Dd = Nd/N the sector densities (so that Da + Dd = 1). The spin overlap 
Q is then just 

X 

For /3</3c = l,Q^OasAf— >cx3 (because the EA order parameter is zero in the paramagnetic phase) while for 
/3 > /9c = 1, (Q^) > as oo, where (•) denotes an average over couplings. So it must be that Da = Dd for /3 < 1 
while Da ^ DdioT P > 1. 

The arguments for the single replica case can be repeated separately within the agree and disagree sectors. Letting 
P± denote the conditional probabilities that Jxoyo = ±N~^^'^ given the other Jxy^s and all ax^s and r^j's, we have 
within either of the two sectors that 

^-ap/VN < e-^'^/^lZ^/Z+p < P+/P_ < e4^/^|Z_/Z+|2 < e^'^/^ (30) 

so that the conditional probability within a single sector Pds for Xoj/o to be doubly satisfied is (1/2) + 0{[3N^^/'^). 
For (3 < Pc, we have Da = 1/2, Dd — 1/2 (in the limit N oo) and so in either sector, double FK percolation is 
approximately a random graph model with A^/2 sites and bond occupation probability (1/2)4/3^ A^~^ ~ P'^{N/2)~^; 
thus double FK giant clusters do not occur for < 1 in the TRFK representation. 

In the CMR representation, blue percolation corresponds to bond occupation probability (1/2)4/3 A^^^/^ = 
f3N^/'^(N/2)~^ and so the threshold for blue percolation is given by /3cmr,p = N^^^^. But now there are two 
giant clusters, one in each of the two sectors, and as noted above they satisfy Da = Dd when /3 < /3c = 1- For fi > (3c, 
Da Dd- Since 

/3iVi/2 ^ oo for /3 > /3c (indeed for any fixed /3 > 0), it follows from random graph theory that each 
giant cluster occupies the entire sector so that Da and Dd are also the cluster percolation densities of the two giant 
clusters. 

In the case of two-replica FK percolation for (3 > (3c, let us denote by Dmax and Dmin the larger and smaller of Da 
and Dd, so that D„iax + Dmin = 1 and D^ax ~ D^ain = Q- Then for (^ > (3c, the bond occupation probability in the 
larger sector is (3'^{N/2)~^ = 2(3"^ Dmax{DmaxN)~^ with 2(3^Dmax > 1 and there is a (single) giant cluster in that larger 
sector. There will be another giant cluster (of lower density) in the smaller sector providing 2/3^ 15mm = /3^(l~Q) > 1- 
Since Q < (/ba, for this to be the case it suffices if for (3 > (3c, 



qEA<l-^. (31) 



The estimated behavior of qsA both as /3 — > H- and as /3 ^ oo [ij suggests that this is always valid. In any case, we 
have rigorously proved that there is a unique maximal density double FK cluster for (3 > (3c- 

An important feature of spin glass order is ultrametricity, which is believed to be true at least for the SK model [s^ , 
although it has not yet been proved rigorously. For the spin overlaps coming from three replicas put into rank 
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order, > Q(2) > Q(3), ultrametricity is the property that Q(2) = Q(3)- The issue of percolation signatures for 
ultrametricity in the SK model is treated at length in |22l |. We only mention here one of those signatures, which 
concerns the four percolating clusters that arise in a CMR representation of three replicas when considering bonds 
that are simultaneously blue both for replicas one and two as well as for replicas one and three. Denoting the four 
cluster densities in rank order as a;(i) > X(2) ^ ^i^i) ^ the percolation version of the ultramctric property is that 
>_X{2) and X(3) = a;(4). The case X{2) = (resp., X(2) > 2:(3)) corresponds to Q(i) > (5(2) (resp, = Q(2))- 
See [22| for more details. 



IV. NUMERICAL METHODS AND RESULTS 



We have carried out numerical simulations to explore the properties of two-replica FK representations described 
in the previous sections for the case of the three-dimensional ± J EA spin glass. The Monte Carlo algorithm that 
we use also makes use of the CMR representation in addition toparallel tempering and Metropolis sweeps. Similar 
methods have been previously applied by Swendsen and Wang [o, [IS 113] and others [sl, [H, . The algorithm is 
described in more detail in [2l|, here we provide some additional details about the CMR component of the algorithm. 
The CMR cluster algorithm alternates between "bond moves" and "spin moves." The bond move takes a pair of 
EA spin configurations {cx} and {tx}, both in the same realization of disorder, {Jxy}, and populates the bonds 
of the lattice with red and blue CMR bonds. Each bond that is satisfied in both replicas is occupied with a blue 
bond with probability T^biue and each bond that is satisfied in only one replica is occupied with a red bond with 
probability Vrcd, defined in (fT5|) and (fT6| . respectively. The bond move is described by the conditional probability, 
T{{ujxy}, {VxyjlWx}, {tx}; {Jxy}) for bond configuration {ujxy}, {Vxy} given spin configuration {ax}, {r^.}. To simplify 
the notation, we omit the lattice indices and the dependence on Jxy in the following equations except where explicitly 
needed. For the bond move we have. 



r(M,M|{a},{r}) = 



i3biuc({c^})i3red({^})A({a}, {t}, {uj})r{{a}, {r}, {77}) 
E{.'}{v} Bbiuc({c.'})i?.cd({^7'})A({a}, {r}, {^'})r({a}, {r}, {r/}) ' 



(32) 



The spin move takes a blue and red CMR bond configuration and produces a pair of spin configurations {(Jx} and {tx}. 
All spin configurations obeying the constraints that blue bonds are doubly satisfied and red bond are singly satisfied 
have the same probability. The spin move is described by the conditional probability, T({ax}, {Tx}\{wxy}, {rixy})- 



^({^},MIM,M) 



A(M,{r},M)r(M,{r},M) 
Ew'},ir'} A({^'}, {r'}, {io})T{{a'}, {r'}, {rj}) ' 



(33) 



A bond move followed by a spin move constitututes a single sweep of the CMR cluster algorithm. We now demon- 
strate the validity of the CMR algorithm by showing that detailed balance and ergodicity are satisfied. Ergodicity 
is clearly satisfied since there is a non-vanishing probability that no red or blue bonds will be created during a bond 
move so that the spin move can then produce any spin configuration. Detailed balance with respect to the equilibrium 
spin distribution can be stated as 



J2 raa}, {r}|K}, {r/'})r({^"}, {rj"}\{a'}, {r'}) x exp 

{uj"}^{r,"} 



pY^JxyKa'y+T'^r;^) 

{xy} 



J2 m^r'}, {r'lIK}, {v"}m{u;"}, {rf}\{a}, {r}) x exp 
{^"}.{v"} 



{^y} 



(34) 



A({a},{r},K'})r(M,{r},{V'}) 



This equation must hold for all choices of {(Jx}, {tx}, {^x} and {r^}. Note that by (|12p and the denominator in 
([5^ is simply the Boltzmann weight for two independent EA spin glasses, cancelling the same factors in the numerator 
of the detailed balance equation. Thus, the LHS of p4|) can be written as 

{.'Wl A({^"}, {r"}, K'})r({a"}, {r"}, {<}) '~ 

Sbiuc({c^})Sred({^})A({a'}, {r'}, K})r({a'}, {r'}, {i^"}) (35) 

Note that this expression is symmetric under the exchange of {ux}, {tx} and {cr^}, {t^} and thus equal to the RHS 
of ([M)) . demonstrating that detailed balance holds. The bond configurations observed after the bond move are, in 



9 



2.0 



1.5 
ul.O 
0.5 



0.0 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

FIG. 1: Average number of wrapping CMR clusters, mJcmr vs. (3 for the 3D EA model. 



the same fashion, easily shown to be equilibrium CMR bond configurations as described by the bond marginal of the 
Edwards-Sokal weight (fT^ . 

Although the CMR algorithm correctly samples equilibrium configurations of EA spin glasses, in 3D it is not 
efficient. To obtain a more efficient algorithm, we also make use of parallel tempering and Metropolis sweeps. The 
CMR algorithm employs two replicas at a single temperature while parallel tempering exchanges replicas at different 
temperatures. Here we use 20 inverse temperatures equally spaced between (3 = 0.16 to (3 — 0.92. The phase transition 
temperature of the system was recently measured as /3c — 0.89 ±0.03 [34|. A single sweep of the full algorithm consists 
of a CMR cluster sweep for each pair of replicas at each temperature, a parallel tempering exchange between replicas 
at each pair of neighboring temperatures and a Metropolis sweep for every replica. 

We simulated the three-dimensional ± J Edwards- Anderson model on skew periodic cubic lattices for system sizes 
6'^, 8'^, 10'^ and 12^. For each size we simulated 100 realizations of disorder for 50,000 Monte Carlo sweeps of which 
the first 1/4 of the sweeps were for equilibration and the remaining 3/4 for data collection. The quantities that we 
measure are the fraction of sites in the largest blue cluster, Ci and second largest blue cluster, C2 and the number of 
blue CMR wrapping cluster, wcmr, and the number of TRFK "wrapping" clusters, wtrfk- A cluster is said to wrap 
if it is connected around the system in any of the three directions. 

Figure [1] shows the average number wcmr of CMR blue wrapping clusters as a function of inverse temperature f3. 
The curves are ordered by system size with largest size on the bottom for the small /3 and on top for large (3. The 
data suggests that there is a percolation transition at some PcMR.p- For (3 > /3cmr,p there are two wrapping clusters 
while for /? < PcMR.p there are none. Near and above the spin glass transition at Pc ~ 0.89 the expected number of 
wrapping clusters falls off but the fall-off diminishes as system size increases. This figure suggests that in the large 
size limit there are exactly two spanning clusters near the spin glass transition both above and below the transition 
temperature. 

Figure [2] shows the fraction of sites in the largest CMR blue cluster, Ci, second largest CMR blue cluster, C2 and 
the sum of the two, Ci + €2- The middle set of four curves is Ci for sizes 6^, 8^, 10^ and 12^^, ordered from top to 
the bottom at /3 = 0.5. The bottom set of curves is C2 with systems sizes ordered from smallest on bottom to largest 
on top at f3 — 0.5. The difference between the fraction of sites in the two largest clusters, Ci — C2 is approximately 
the spin glass order parameter. As the system size increases, this difference decreases below the transition suggesting 
that Ci = C2 for /3 < /?c in the thermodynamic limit. On the other hand, the sum of the two largest clusters is 
quite constant independent of system size. Near the transition, approximately 96% of the sites are in the two largest 
clusters. 

The large fraction of sites in the two largest clusters makes the CMR cluster moves inefficient. If all sites were in 
the two largest clusters then the cluster moves would serve only to flip all spins in one or both clusters or exchange 
the identity of the two replicas. Equilibration depends on the small fraction of spins that are not part of the two 
largest clusters. 

Figure [3] show the average number of wrapping TRFK clusters wtrfk as a function of inverse temperature. The 
largest system size is on the bottom for the small {3 and on top for the large /3. As for the case of CMR clusters, 
the data suggests a transition at some /3trfk,p from zero to two wrapping TRFK clusters. Although the number of 
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FIG. 2: Ci (middle set), C2 (bottom set) and Ci +C2 (top set) vs. /3 for the OMR graphical representation for the 3D EA model. 
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FIG. 3: Average number of doubly occupied wrapping Fortuin-Kasteleyn clusters, wtrfk vs. /3 for the 3D EA model. 

TRFK wrapping clusters is significantly less than two for all (3 and all system sizes, the trend in system size suggests 
that it might approach two for large systems and /3 > Ptrfk.p- 

The percolation signature for both CMR and TRFK clusters is qualitatively similar in three dimensions. In both 
cases two giant clusters with opposite values of the local order parameter appear at a temperature substantially above 
the phase transition temperature. In the high temperature phase, the two giant clusters have the same density and 
the phase transition is marked by the onset of different densities of the two clusters. The numerical evidence, however, 
suggests that the transition from equal to unequal giant cluster densities is quite broad for the small system sizes 
explored here. 



DISCUSSION 



We begin by noting that our results, even for the SK infinite-range spin glass, are not used to prove a phase 
transition (which has already been proved using other techniques 0). Rather, we use as an input the fact that such 
a transition exists and that the EA order parameter is zero above Tc and nonzero below, and then show that the 
transition coincides exactly with the onset of a density difference in doubly occupied FK clusters. Numerical evidence 
indicates that something similar is occurring in the EA model. Together these complementary approaches shed light 
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on the nature of the spin glass transition, especially from a geometric viewpoint, and suggest a possible framework 
through which a spin glass phase transition in realistic spin glass models can be finally proved. 

We now summarize our results. We have introduced a random cluster approach for studying phase transitions and 
broken symmetry in spin glasses, both short- and infinite-range. We have shown that, unlike for ferromagnetic models, 
single FK percolation is a necessary but not sufficient condition for broken spin fiip symmetry. However, double FK 
percolation (with a unique largest cluster) is sufficient and probably necessary. (More precisely, it is necessary in the 
SK model, because otherwise the EA order parameter is zero; and it is probably necessary in the EA model.) 

In the SK model, there is a difference above Tc between the CMR and TRFK approaches, but not below. In the 
former, there already is double percolation below /3c, above an onset inverse temperature of /3cmr,p = N~^/^, below 
which there are no giant clusters. Between this inverse temperature and Pc (which equals one in the class of models 
we study) there are exactly two giant clusters of equal density, which become 1/2 in the N ^ oo limit. 

More importantly, there is a second transition occurring at exactly the SK spin glass critical value /?c = 1. Above 
this threshold, the two giant clusters take on unequal densities, whose sum is one (i.e., every bond and spin belong to 
one of the two giant clusters). It could be the case at some even higher /3 there is only a single giant cluster, but our 
methods so far are unable to determine whether this is the case. 

For TRFK percolation there are no giant clusters above /3trfk,p = 1- For /3 > 1 there are exactly two giant clusters 
with unequal densities, and the picture then becomes similar to that of the CMR representation. 

The numerical simulations of the 3D EA model suggest a scenario similar to what we find in the SK model. We 
observe a sharp percolation transition for both the CMR and TRFK representations at a temperature well above the 
spin glass transition temperature. At this percolation transition, two giant clusters form with opposite values of the 
EA order parameter. These clusters are nearly equal in density and together occupy most of the system. As the 
temperature is decreased toward the presumed location of the spin glass transition, the density of the two clusters 
becomes increasingly unequal. Although this transition in the density of the two largest clusters appears quite rounded 
for the small systems investigated numerically, we believe it is sharp in the thermodynamic limit. The results for both 
the SK model and the connection between the EA order parameter and the density diflference between giant clusters 
strongly suggests that the spin glass transition in finite dimensions is marked by the onset of this density difference 
in both the CMR and TRFK representations. These results provide an interesting geometric avenue for investigating 
the phase transition in the EA model and related models with frustration such as the random bond Ising model or 
the Potts spin glass. 
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